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We present an extensive analysis of gauge-mediated supersymmetry breaking models with minimal 
and non-minimal flavour violation. We first demonstrate that low-energy, precision electroweak, and 
cosmological constraints exclude large "collider- friendly" regions of the minimal parameter space. 
We then discuss various possibilities how flavour violation, although naturally suppressed, may 
still occur in gauge-mediation models. The introduction of non-minimal flavour violation at the 
electroweak scale is shown to relax the stringent experimental constraints, so that benchmark points, 
that are also cosmologically viable, can be defined and their phenomenology, i.e. squark and gaugino 
■ production cross sections with flavour violation, at the LHC can be studied. 
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43 ' I. INTRODUCTION 

£LT Weak scale supersymmetry (SUSY) remains a both theoretically and phenomenologically attractive extension of 
the Standard Model (SM) of particle physics [H,[2]- Apart from linking bosons with fermions and unifying internal and 
external (space-time) symmetries, SUSY allows for a stabilization of the gap between the Planck and the electroweak 
scale and for gauge coupling unification at high energies. It appears naturally in string theories, includes gravity, and 
contains a stable lightest SUSY particle (LSP) as a dark matter candidate. Spin partners of the SM particles have not 
yet been observed, and in order to remain a viable solution to the hierarchy problem, SUSY must be broken at low 
energy via soft mass terms in the Lagrangian. As a consequence, the SUSY particles must be massive in comparison 
to their SM counterparts, and the Tevatron and the LHC will perform a conclusive search covering a wide range of 
masses up to the TeV scale. After the discovery of SUSY particles, the revelation of the underlying SUSY-breaking 
• ' mechanism will be one of the key challenges in the experimental high-energy physics program. 

In gauge-mediated supersymmetry-breaking (GMSB) models, SUSY is broken in a secluded sector at a scale (F) 
related to the gravitino mass by m,Q — (F) / (y/3Mp) , where Mp is the reduced Planck mass. The breaking is 
mediated to the visible sector of squarks, sleptons, gauginos and gluinos through a gauge-singlet chiral superfield 
S and n q quark- like and ni lepton-like messenger fields [3, 0, B @] • The superfield S is characterized by its scalar 
,£h ! an d auxiliary components, which overlap with the gravitino and acquire vacuum expectation values (S) and (Fs), 
respectively. Yukawa couplings of the messengers to the superfield S then induce masses of order M mcs ~ (S) for the 
messengers. Gauginos and sfermions acquire masses through ordinary gauge interactions with messengers through 
one- and two-loop self-energy diagrams, respectively. In these scenarios, the lightest SUSY particle is always the 
gravitino, which is thus a natural candidate for the dark matter in our Universe. Besides M mcs , n q , and ni, minimal 
GMSB scenarios are determined by the ratio of the two Higgs vacuum expectation values, tan/3, the sign of the off- 
diagonal Higgs mass-parameter /i, and by the auxiliary vacuum expectation value (Fs). The latter is related to the 
mass splitting of the messenger fields and is considerably smaller than both the squared mass scale of the messenger 
fields, (S) 2 , and the fundamental SUSY-breaking scale, (F). It is usually re-expressed in terms of an effective SUSY- 
breaking scale, A = (Fs)/(S). An additional free parameter is the gravitino mass, mg, which is, however, constrained 
by the fact that the gravitino relic density fi^/i 2 has to agree with the current WMAP limits and that the abundances 
of the light elements should be correctly described, i.e. the next-to-lightest SUSY particle (NLSP) must not decay too 
quickly. 

GMSB is an attractive scenario regarding the so-called SUSY flavour problem. SUSY is usually broken within a 
few orders of magnitude of the weak scale, whereas the unrelated flavour-breaking scale can be chosen much higher. 
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This avoids important flavour-violating terms in the SUSY-breaking Lagrangian and leads to approximately flavour- 
conserving mass matrices at the low-energy scale and good agreement with measurements of flavour-changing neutral 
current observables. However, several possibilities reintroducing flavour-violating terms in the Minimal Supersym- 
metric Standard Model (MSSM) with GMSB have been pointed out @,0,[1]- 

For example, mixing between messenger 
and matter fields may lead to important flavour violations in the squark and slepton sectors. 

In SUSY models with non- minimal flavour violation (NMFV), the flavour- violating off-diagonal terms Aij of the 
squared sfermion mass matrices, where i,j — L,R refer to the helicities of the (SM partners of the) sfermions, are 
conveniently considered as arbitrary parameters. Stringent experimental constraints are then imposed by precise 
measurements of K° — K° and B° — B° mixing, the first evidence of D° — D° mixing, and rare decays dUUH. The 
minimal GMSB model obviously relies on constrained minimal flavour violation (cMFV), where all the flavour- violating 
elements A^ are neglected. Recently, possible effects of non-minimal flavour violation on the experimentally allowed 
minimal supergravity (mSUGRA) parameter space have been investigated, and all squark and gaugino production 
cross sections and decay widths have been recalculated including both helicity and flavour mixing in the squark sector 
fl2l ]. The aim of this work is to extend this study to GMSB scenarios, to evaluate the experimental constraints, discuss 
the role of flavour violation, and make numerical predictions for squark- and gaugino-production cross sections at the 
LHC. 

This paper is organized as follows: In Sec. UH we impose the current experimental constraints on the minimal 
GMSB models. We show that these scenarios are strongly disfavoured due to the very stringent constraint coming 
from the rare b — > decay. However, the latter can be relaxed by introducing NMFV in the squark sector, as shown 
in Sec. IIII1 allowing us to define benchmark points for NMFV GMSB scenarios. Sec. [IV] is devoted to the discussion 
of cosmological implications on the gravitino mass in our scenarios. In Sec. El we present numerical predictions for 
squark and gaugino hadroproduction cross sections at the LHC. Our conclusions are given in Sec. IVII 



II. EXPERIMENTAL CONSTRAINTS ON GMSB MODELS WITH MINIMAL FLAVOUR VIOLATION 

In the absence of experimental evidence for Supersymmetry, a large variety of data can be used to constrain the 
parameter space of the MSSM. Sparticle mass limits can be obtained from searches of charginos (m~± > 150 GeV 
from DO), neutralinos (m^o > 93 GeV in GMSB from the combination of LEP2 results), gluinos (rrig > 195 GeV from 
CDF), stops {m$ > 95 ... 96 GeV for neutral- or charged-current decays from the combination of LEP2 results), other 
squarks (m§ > 300 GeV for gluinos of equal mass from CDF), and gravitinos (mg > 1.3 • 10 -5 eV for rriq = m g = 200 
GeV) at colliders [H,[l3. 

Cosmological, electroweak precision, and low energy observables can be used to put additional constraints on the 
SUSY parameter space. The theoretically robust inclusive branching ratio 

BR(6^s 7 ) = (3.55 ± 0.26) ■ 10" 4 , (1) 

obtained from the combined measurements of BaBar, Belle, and CLE O fla| , can be confronted to theoretical pre- 
dictions including two-loop QCD and one-loop SUSY contributions (T5 , Il7||. Squarks contribute here already at the 
one-loop level, as do the SM contributions. A second observable, sensitive to the squark-mass splitting within isospin 
doublets, is the electroweak p-parameter with 

E z (0) 3V(0) 

Ap = § — ' ( 2 ) 

where mz,w an d £z,w(0) denote the Z- and VF-boson masses and self-energies at zero momentum, respectively. New 
physics contributions are constrained by the latest combined electroweak precision measurements to T = —0.13 ±0.11 
or 

Ap = -aT = (1.02 ± 0.86) • 10" 3 (3) 

for a(mz) = 1/127.918 This value is compared to theoretical calculations including SUSY two-loop corrections 
fl8l ]. A third variable sensitive to new physics loop contributions is the anomalous magnetic moment of the muon, 
for which we require the SUSY contribution a^ USY , known up to the two-loop level [H, 0]> to c l° se the gap between 
recent BNL experimental data and the SM prediction [l3j . 

Aa p = (29.2 ±8.6) • 10~ 10 . (4) 

Note that the one-loop SUSY contributions are approximatively given by [2l[ 

a susY,i-ioo P ^ 13. 10- W (^^V tan/3 sgnOi), (5) 
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FIG. 1: The A-M mos planes for /i > and different values of tan/3 and AT mes , assuming constrained minimal flavour violation. 
Dark (blue) regions are excluded by the constraint coming from the b — > s'y branching ratio. The regions where A > M mes does 
not allow for physical solutions of the RGEs. 



if the relevant SUSY particles have masses of the order of Msusy- As a consequence, negative values of /i then 
increase, not decrease, the disagreement between the experimental measurements and the theoretical value of a^, so 
that the region /i < is strongly disfavoured in all SUSY models. In addition, this region is also virtually excluded 
by the b — > sj constraint at the 2tr confidence level. We therefore restrict ourselves to positive values of \i throughout 
this analysis. 

The above experimental limits are imposed at the 2er confidence level on the minimal GMSB model with /i > 
and four free parameters A, M mes , 7V mes = n q — ni, and tan/?. The renormalization group equations (RGEs) are 
solved numerically to two- loop order using the computer programme SPheno 2.2.3 [22j , which computes the soft 
SUSY-breaking masses at the electroweak scale with the complete one-loop formulas, supplemented by two-loop 
contributions for the neutral Higgs bosons and the /i-parameter. We then diagonalize the mass matrices and compute 
the electroweak precision and low-energy observables with the computer programme FeynHiggs 2.6.4 [23| . For the 
SM input parameters, i.e. the masses and widths of the electroweak gauge bosons and quarks, the ang les of the 
CKM-matrix and its CP- violating phase, and the Fermi coupling constant, we refer the reader to Ref. 

In Fig.[TJ we show typical scans of the minimal GMSB parameter space in A and M mes for different values of tan p 
(15, 30, and 50) and N mcs (1 and 3). The six panels reveal that these scenarios are stro ngly disfavoured by the 
measurements of the b — > branching ratio. In particular, the Snowmass benchmark points [24| SPS 7 (A = 40 TeV, 
M mcs = 80 TeV, tan/3 = 15, n > 0, and N mes = 3) and SPS 8 (A = 100 TeV, M mcs = 200 TeV, tan/3 = 15, n > 0, 
and Ames = 1) lead to values of BR(b — > s-f) — 6.97 • 10~ 4 and 6.77 • 10 -4 , which are both excluded beyond the 5a 
level, even if both of these points lie well within 2er of the experimentally allowed range for the anomalous magnetic 
moment of the muon with a^ USY = 22.8 • 10~ 10 and a^ USY = 16.31 • 10~ 10 . Note that the regions with A > Af mcs are 
theoretically excluded, since they do not allow for physical solutions of the RGEs. 

Recently, a detailed study of electroweak precision observables, including scenarios with minimal GMSB, has been 
performed ^25J. Scanning also over tan/3 and allowing for higher values of N mcs < 8, the authors show that experi- 
mentally favoured scenarios can be achieved at low messenger scales, which, however, implies a certain amount of fine 
tuning at the weak scale, e.g. in the Higgs sector. Note that for N mcs > 8 problems with perturbativity of the gauge 
interactions arise at very high scales 
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III. GMSB MODELS WITH NON-MINIMAL FLAVOUR VIOLATION 



The minimal GMSB is known to suppress flavour-changing neutral currents as suggested by measurements and thus 
to avoid the SUSY "flavour problem" , which arises naturally in models where SUSY-breaking is mediated by gravity. 
Models beyond the minimal GMSB can, however, reintroduce flavour-breaking terms at the electroweak scale. In 
this Section, we first review flavour violation in the MSSM, present its implementation at the electroweak scale, and 
elaborate on different non-minimal GMSB models including flavour violation. We then re-analyze the parameter space 
and show how NMFV can provide a way to relax the stringent constraints challenging the minimal GMSB models. 



A. Theoretical Framework 



In constrained minimal flavour violation (cMFV) SUSY models, the only source of flavour violation arises through 
the rotation of the quark interaction eigenstates into the basis of physical mass eigenstates, where the Yukawa matrices 
are diagonal, as in the SM, and the flavour-violating entries of the squark mass matrices are neglected both at the 
SUSY-breaking and the weak scale. In SUSY with NMFV, these flavour-violating entries A?? are considered as free 
parameters. The squared squark mass matrices are then given by 



/ Ml AH A" 
Afl Ml Af L 
Afl Afl Ml 


mxXx A\l \ 
AJ£ m 2 X 2 Af R 
A%Z Afl m 3 X 3 


mi Xi A R \ 
Af^ m 2 Y| Af L 
V A i« A|^ m 3 Y 3 * 


M 2 Rl AH &RR 
Ajjji Ml A» 

A nn M k J 



M? 

q 



where M\ and M R denote their usual diagonal entries, 

M 2 Lk = M% k + ml + cos 2/3 m| (T%-e k sm 2 w ) 
M R k = M u k ,D k + m t + cos 2/3 m z e k sin 2 W) 
and hclicity mixing is generated by the elements 



(6) 



(7) 

(8) 



X q = 



At-n 




for up — type squarks, 
for down — type squarks. 



(9) 



Here, 9w is the electroweak mixing angle, Mq^d are the usual SUSY-breaking squark masses, A q is the trilinear 
coupling, and mf., T^, and denote the mass, weak isospin, and electric charge of the quark the index k referring 

to the (s)quark generation. The flavour- violating elements A|? are usually normalized to the diagonal entries [Ioj |. 



A" = 



Xf, Mi Mi , . 

11 L q Jq' ■ 



(10) 



so that NMFV is governed by 24 arbitrary complex dimensionless parameters Xfj . 

The diagonalization of the mass matrices M| and Mj requires the introduction of two additional 6x6 matrices 
R u and R d , relating the hclicity and flavour eigenstates to the physical mass eigenstates through 



R u {u L ,c L ,t L ,u R , c R , t R ) T , 



({h,«2,U3,U4,U 5 ,M 6 ) T 

(di,J 2 ,rf3,rf4,d5,rf6) T = R d (dL 1 s L ,b L ,d R ,s R ,b R ) T 



(11) 

(12) 



By convention, the squark mass eigenstates are labeled according to < ... < for q = u,d. For a detailed review 
of flavour violation in the MSSM see e.g. Ref. [26|. We stress that we do not employ the mass-insertion approximation 

with its perturbative expansion in the parameters A?J , but rather perform the diagonalization of the squark-mass 
matrices numerically. Relatively strong constraints on NMFV SUSY models can be obtained from low-energy and 
electroweak precision observables, e.g. upper limits from the neutral kaon sector, on B- and D-meson oscillations, 
various rare decays, and electric dipole moments. In several publications [ToL ITD. \2l\ . rather complete analyses have 
been presented, pointing out that the down-squark sector is particularly constrained from K- and _B-physics processes 
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with external down-type quarks and that within the mass-insertion approximation the only substantial mixing in the 
squark sector occurs between the second and third generations in the left-left and right-right chiral sectors. Note that 
the latter is suppressed in gravity mediation models by the scaling of the corresponding entries A?? with the SUSY 
breaking scale, while in gauge mediation models the mixing in the left-right chiral sector, induced by A-terms, is 
small. The up-squark sector is in general less experimentally constrained. This situation may change once additional 
information from neutral and charged Higgs production and decay becomes available [l6l . [28j . In our analysis, we 
apply SU(2) gauge invariance to the left-chiral sector and take implicitly into account the above mentioned constraints 
by restricting ourselves, also for the sake of simplicity, to the case of two real NMFV parameters, 

A LL = Aft ~ Aft < 0.2 and A RR = A&. ^ < 0.2, (13) 
while all other A? 9 are zero. 

3 

Although the gauge interactions are flavour-blind, it has been shown that there are several possibilities for flavour 
violation in both the squark and slepton sectors to arise within GMSB models. For example, for very high messenger 
scales M mcs > 10 15 GeV gravity is no longer negligible with respect to gauge interactions @. As a consequence, 
flavour-violating terms are reintroduced through gravity mediation as in mSUGRA models. However, scenarios with 
such high messenger scales are rather unattractive from a phenomenological point of view due to the resulting very 
high SUSY masses. Second, flavour violation can also be induced from heavy right-handed neutrinos participating in 
leptogenesis [7[ • If these are lighter than the messenger scale, flavour off-diagonal mass terms are introduced into the 
slepton mass matrices. A third possibility might be to consider broken messenger number invariance, that implies 
that the lightest messenger is not stable and introduces flavour-violating terms in the Lagrangian at the weak scale 
Q . A disadvantage of this model is that the now unstable lightest messenger may not be a viable candidate for cold 
dark matter in the case of a very light gravitino, which cannot account for the observed relic abundance. 

For our study, we focus on the model proposed in Ref. Q , based on the introduction of a mixing between messenger 
and matter fields. In the case of fundamental messenger multiplets belonging to 5 and 5 representations of SU(5), the 
messengers carry quantum numbers of left-handed leptons and right-handed down-type quarks. As a consequence, 
flavour violation is introduced into the chiral sectors of right-handed sleptons and left-handed (up- and down-type) 
squarks. In our analysis of squark flavour violation, this corresponds to including a variation of the parameter All, 
while Arr is set to zero. An alternative model with antisymmetric messenger multiplets belonging to 10 and 10 
representations can also be considered. In this case, the messengers share quantum numbers with right-handed 
leptons, left-handed up- and down-type quarks, and right-handed up-type quarks, leading to flavour mixing for left- 
handed sleptons as well as for both left- and right-handed up- and down- type squarks. Note that, in this antisymmetric 
scenario, flavour mixing in the sector of right-handed down-type squarks may be parameterized independently of the 
other squarks. In our analysis, however, we use for simplicity the same flavour violation parameter All = Arr for 
both chiral sectors. We should stress that in both scenarios flavour violation is completely governed by the parameter 
All- 



B. Scans of the Parameter Space and Benchmark Points 

We now re-investigate the constraints discussed in Sec. [TT] and include non-minimal flavour violation as discussed 
above. Allowed regions for the parameters A, M mcs , A mos , tan/3, sgn(/i), All, and Arr are obtained by explicitly 
imposing the constraints from b — > sj, Ap, and a^, that are sensitive to flavour-violating terms. In particular, 
the branching ratio BR(6 — > s"f) is directly affected by the allowed squark mixing between the second and third 
generation. Squarks enter the calculation at the one-loop level, as do the SM contributions, so that the dependence 
on the NMFV-parameters is rather important. A second important consequence of NMFV in the MSSM is the 
generation of a large splitting between squark-mass eigenvalues, which directly influences the electroweak precision 
variable Ap. Concerning a^ USY , the squark contributions are suppressed with respect to the slepton contributions, 
so that its dependence on flavour violation is less important. The renormalization group running is again performed 
with SPheno 2.2.3. The flavour violating terms All and Arr are included in the squark sector at the weak scale as 
discussed in Sec. lIII Al before diagonalizing the mass matrices and computing the low-energy and electroweak precision 
observables with FeynHiggs 2.6.4. 

In Figs. [5]- [7| we show scans of the A-M mes plane for p > and the same values N mes — 1, 3 and tan/3 = 15, 30, 50 
as in Sec.|lIJ The region favoured by the anomalous magnetic moment of the muon (light/grey) is quite insensitive 
to variations of the parameter All, since dominant SUSY effects come from induced quantum loops of gauginos 
and sleptons, while squarks contribute only at the two-loop level, which reduces the dependence on squark flavour 
violation. As expected, the b —* sj excluded region (dark/blue) depends strongly on flavour mixing, while the 
constraint coming from Ap does not play a role for the moderate SUSY masses relevant to our region of interest, so 
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FIG. 2: The A-M mes planes for Y mcs = 1, tan/3 = 15, fi > 0, and All = 0.1,0.15 and 0.2. We show a M favoured (light/grey) 
and b — > S7 excluded (dark/blue) regions of the GMSB parameter space with non-minimal flavour violation in either the 
left-left chiral (Arr = 0, top) or both the left-left and right-right chiral (Arr = All, bottom) squark sectors. The region where 
A > M mcs does not allow for physical solutions of the RGEs. 



that the corresponding excluded regions are not shown. Note that the difference between the two flavour violation 
scenarios considered here is relatively small, as can be seen in Figs.[2J-[7] It becomes clear that, if we allow for flavour 
mixing between the second and third generation squarks, windows in the parameter space both favoured by a M and 
not excluded by the stringent constraint from b — > sj make their appearance for small and moderate SUSY masses. 

Within these regions, we propose six benchmark scenarios permitting non- minimal flavour violation and not yielding 
too high SUSY masses ("collider-friendly"), so that possible SUSY signals should be observable at present and/or 
future hadron colliders. Our choices are presented in Tab.[H labeled starting at the point E due to our four benchmark 
proposals for mSUGRA scenarios including NMFV Note that, in contrast to the mSUGRA case, these scenarios 
are not valid assuming cMFV (All = Arr = 0), so that we indicate the allowed ranges for our flavour mixing 
parameter All- 

Starting with tan/3 = 15 and N mcs = 1 (see Fig. 12]), we choose our benchmark point E in the region both favoured 
by the electroweak precision constraints and corresponding to rather light SUSY particles. As for any GMSB scenario, 
the gravitino is the lightest SUSY particle (LSP). The next-to-lightest SUSY particle (NLSP) is the lightest neutralino 
with m^o = 95.4 GeV, but the three lightest charged sleptons are very close with similar masses around 100 GeV. 
The other sleptons, sneutrinos, and gauginos have moderate masses of about 150 - 300 GeV, while the squarks and 
gluino are quite heavy with masses lying in the range of 700 - 800 GeV. However, they are much lighter than those 
corresponding to the point SPS 8 with its larger values of A and M mes , lie well above the experimental limits obtained 
from direct searches assuming cMFV, and are experimentally accessible at the LHC. 

The point F (see Fig. [3]) differs very little from the point SPS 7, with the SUSY-breaking scale A shifted from 40 
to 30 TeV, so that it now lies in the preferred region with respect to the b — ► sj constraint. As for SPS 7, the three 
lightest sleptons have masses around 100 GeV, the lightest being the stau with = 90.7 GeV. The other sleptons, 
sneutrinos, and gauginos are a bit heavier (120-200 GeV), and the squarks and gluino are rather heavy (600-700 GeV). 

The points G, H, I, and J (see Figs. 2][7]) all have a stau NLSP with a mass between 99 and 160 GeV. The main 
difference in their spectra is the number and nature of the particles that are closest in mass to the NLSP. For the 
point G, these are two sleptons and the lightest neutralino, whereas for the point H these are only the two sleptons. 
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N„=3, tar$=15,)i>0,\ =0,\ LL =0.1 N ms =3, tar#=15,\i>0,\ =0,l LL =0.15 N m =3, tan(,= 15,\i>0,\ =0,\ LL =0.2 
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FIG. 3: Same as Fig. [2]for N mes = 3 and tan (3 = 15. 
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FIG. 4: Same as Fig.[2]for iV mcs = 1 and tan/3 = 30. 
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N„=3,tar$=30,)i>0,\ =0,\ LL =0.1 N mes =3, tar$=30,\i>0,X =0,l LL =0.15 N ms =3, tar(,=30,\i>0,\ =0,\ LL =0.2 




50 100 150 200 50 100 150 200 50 100 150 200 

M„,„ (TeV) M™, (TeV) M m „ (TeV) 



FIG. 5: Same as Fig. [5] for iV mcs = 3 and tan /3 = 30. 
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FIG. 6: Same as Fig. H for iY mcs = 1 and tan/3 = 50. 
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N mas =3,tar$=50,)i>0,\ =0,\ LL =0.1 N mes =3, tar$=50,\i>0,X =0,l LL =0.15 N ms =3, tarf,=50,\i>0,l. =0,\ LL =0.2 
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FIG. 7: Same as Fig. [5] for iV mes = 3 and tan /3 = 50. 



TABLE I: GMSB benchmark points allowing for non-minimal flavour violation in the left-left (Arr = 0) or both the left-left 
and right-right (Arr = All) chiral squark sectors. We also indicate the allowed range for the NMFV-parameter All, the nature 
of the next-to-lightest SUSY particle (NLSP), and the closest SPS benchmark point (if relevant), which are valid for both 
flavour-violating scenarios. 





A [TeV] M mes [TeV] A mcs tan/3 sgn( Al ) 


All 


NLSP 


SPS 


E 
F 


65 90 1 15 + 
30 80 3 15 + 


[0.14, 0.20] 
[0.12, 0.18] 


X? 
fi 


8 
7 


G 
H 


100 110 1 30 + 
45 100 3 30 + 


[0.14, 0.20] 
[0.12, 0.18] 


h 
h 




I 

J 


130 140 1 50 + 
60 100 3 50 + 


[0.14, 0.20] 
[0.14, 0.20] 


fi 
fi 





The points I and J do not have any particles close to the NLSP in mass. For the four points G, H, I, and J, the other 
sleptons and gauginos are rather light (200-600 GeV), while the squarks and the gluino are very heavy (1-1.5 TeV). 

We now study in detail the dependence of the electroweak precision and low-energy observables as well as the mass 
spectra of the points E, F, G, H, I, and J. In Figs.[S]-[T31 we show the corresponding branching ratio BR(6 — > sj) (top 
left) and the observable Ap (bottom left) as a function of the NMFV-parameter All- We include both flavour violation 
scenarios, the one for fundamental messengers with Arr = as well as the one with antisymmetric messengers, where 
Arr = All- As already mentioned, the leptonic observable a M depends only very weakly (at the two-loop level only) 
on the squarks. As a consequence, we find values of a M independent of All for our six benchmark scenarios, which 
are a^ USY = 37.7 • lO" 10 , 41.3 • lO" 10 , 31.4 • 10" 10 , 36.6 • 10" 10 , 31.8 • lO" 10 , and 34.2 • 10" 10 for the points E, F, G, 
H, I, and J, respectively. These values lie well within 2er of the experimentally favoured range of Eq. ([4]), and even 
within lcr for the points E, G, H, I, and J. For the inclusive branching ratio BR(6 — > sj), the experimentally allowed 
range within 2a is indicated by two horizontal dashed lines. The good agreement between the measurements and 
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FIG. 8: Dependence of the precision variables BR(fe — > sj) and Ap as well as of the up- and down-type squark masses on the 
NMFV-parameter All for flavour mixing in the left-left (Arr = 0) or both the left-left and right-right (Arr = All) chiral 
squark sectors for our benchmark scenario E. The experimentally allowed ranges within 2a are indicated by horizontal dashed 
lines. The vertical dotted lines indicate the allowed range for All with respect to the most stringent constraint from b —* s"f. 
For All = Arr > 0.9 no physical solutions are possible. 



the two-loop SM prediction in combination with the strong dependence of the SUSY contribution on squark flavour 
mixing only leave two allowed narrow intervals for our flavour violation parameter, one being at relatively low values 
of All ~ 0.15, the second one at higher values of All ~ 0.5 — 0.7. It is well known that the latter is disfavoured 
by b — > s(i + (j.~~ data [29| . The remaining one is indicated by vertical dotted lines. Note that the difference between 
the two scenarios is small for the relevant values of All ^ 0.2. Concerning the observable Ap, the difference between 
the two considered flavour mixing scenarios is not visible, so that only one curve is shown. Again, the horizontal line 
indicates the favoured range within 2a, where only the upper limit is visible on our logarithmic scale. In contrast 
to BR(6 — > S7), here the relatively large experimental errors allow for values of All ^ 0.3 — 0.6, depending on the 
benchmark point. The vertical dashed lines indicate the allowed range for All with respect to the more stringent 
constraint coming from b — > 57. 

The difference between the two flavour violation scenarios becomes more obvious when we study the squark mass 
eigenvalues. The up- and down-type squark masses are shown as a function of the flavour mixing parameter All 
in the centre and right upper panels of Figs. [8] -[13] for mixing with fundamental messengers (Arr = 0) and in the 
centre and right lower panels for antisymmetric messengers (Arr = All)- We observe here the same level-reordering 
phenomenon between neighbouring states as already in the case of minimal supergravity, discussed in Ref . [l2j . With 
increasing flavour violation the mass splitting between the lightest and heaviest mass eigenstates becomes larger, while 
the intermediate squark masses are practically unchanged. At the points, where two levels should cross, we observe 
so-called "avoided crossings" of the mass eigenvalues. This phenomenon of level-reordering is due to the fact that 
the mass matrices depend on on single real parameter All- Unfortunately, many "avoided crossings" lie below, but 
some also within the allowed ranges of the flavour- violating parameter All, indicated by green vertical lines for each 
of the six benchmark scenarios. The level-reordering phenomenon is of similar importance for up-type and down-type 
squarks. Concerning the difference between our two implementations of flavour violation in the squark matrices, we 
observe an important splitting for only the lightest and heaviest eigenstates in the case of flavour mixing only in 
the left-left chiral sector. In contrast, for flavour violation in both the left-left and right-right chiral squark sectors, 
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FIG. 9: Same as Fig. [S]for our benchmark scenario F. 




FIG. 10: Same as Fig. [8]for our benchmark scenario G. 
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FIG. 11: Same as Fig. [5]for our benchmark scenario H. 




FIG. 12: Same as Fig. [8]for our benchmark scenario I. 
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the two lightest and two heaviest mass eigenvalues give rise to an important splitting, while only the remaining two 
masses are practically independent of All ■ This is a direct consequence of the fact that wc have introduced additional 
flavour mixing in two distinct sectors of the squark mass matrices, and this will influence the squark and gaugino 
production cross sections presented in Sec. [Vj Note also that in the case of flavour mixing only in the left-left chiral 
sector, "avoided" level crossings occur among the q\.2, 53,4, and q^^ mass eigenstates, whereas in the case of flavour 
mixing in both the left-left and right-right chiral squark sectors, we rather observe the mass flips among the 52,3 and 
94 5 mass eigenstates, respectively. 



IV. COSMOLOGICAL CONSTRAINTS 



To be cosmologically viable, a supersymmetric model should include a convincing candidate for the cold dark matter 
(CDM) in our Universe. This particle has to be stable, electrically neutral, and a colour singlet [13, [M|. Furthermore 
its relic density has to lie within the range 

0.094 < n CDM ft 2 < 0.136 (14) 

at 95% (2a) confidence level. Here, h denotes the present Hubble expansion rate in units of 100 kms -1 Mpc -1 . This 
limit has been obtained from the three-year data of the WMAP satellite, combined with recent SDSS and SNLS 
survey and Baryon Acoustic Oscillation data, and interpreted within an eleven-parameter inflationary model [32| . 
which is more general than the usual six-parameter "vanilla" concordance model of cosmology. Note that this range 
is well compatible with the older, independently obtained range of 0.094 < f^cDM^ 2 < 0.129 of Ref. (33[. 

A natural candidate in SUSY models is the lightest supersymmetric particle (LSP), which is the gravitino in GMSB 
models. Depending on its mass, the gravitino can account either for cold (mg > 100 keV), for warm (1 keV < < 
100 keV), or for hot (ttiq < 1 keV) dark matter. Today's gravitino abundance in the Universe has two contributions. 
First, gravitinos are produced by thermal scattering in the very early Universe. The corresponding energy density 
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FIG. 14: The cosmologically favoured regions in the itiq-Tb. plane for our benchmark scenarios E, F, G, H, I, and J with respect 
to WMAP data (dark/blue), the NLSP lifetime (vertical red line), and leptogenesis (horizontal green line). Also indicated are 
the limits between gravitino warm (WDM) or cold (CDM) dark matter and the regions where the gravitino would not be the 
LSP. 



T R \ /100 GeV\ 



fi£fc a ~ 0.27 — ^ (15) 

Vl0 10 GeV/V m 6 J \1 TeV J V ' 

involves the gluino mass m~ g at low energy and the reheating temperature Tr. The latter is the temperature of the 
Universe after inflation, for which at present no stringent constraints exist. Values of Tr > 10 9 GeV are preferred 
in scenarios that feature leptogenesis in order to explain the cosmic baryon asymmetry 37]. As the resummation 
method leading to Eq. (15) may become unreliable below Tr ~ 10 7 GeV, we use the more accurate result given in 
[HI and do not make practical use of the low-temperature region (see below). Note that the thermal gravitino relic 
density may also be affected by late-time entropy production coming, e.g., from the decay of the messengers. 

Second, there is non-thermal production through decay of the next-to-lightest supersymmetric particle (NLSP) into 
the gravitino. As each NLSP will decay into its Standard Model partner and one gravitino, the resulting gravitino 
energy density can be obtained through 

fi non-t V = _^L_ ^th /l 2 ) (16) 
G "1NLSP 

where ^nlsp^ 2 * s the thermal freeze-out relic density the NLSP would have if it did not decay. Note that for low values 
of m,Q and/or high reheating temperatures Tr thermal production dominates, whereas for high values of tjIq it is 
negligible with respect to the contribution from NLSP decay. The NLSP would-be relic density ^nlsp^ 2 is calculated 
by solving the Boltzmann equation and can be evaluated numerically for any type of NLSP using the public code 
micrOMEGAs [3§], whereas the programme DarkSUSY [4(| is only adapted to the neutralino case. We therefore always 
use micrOMEGAs in this work. 
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Another constraint for scenarios with gravitino dark matter arises from the fact that the NLSP spoils the abundances 



of light elements in our Universe, if it does not decay rapidly enough [35j . The lifetime of a supersymmetric particle 
decaying into its Standard Model partner and a gravitino is given by the inverse of the corresponding decay rate. 
Neglecting here flavour violation and any SUSY particle mixing, which have only little impact [12j , we obtain for the 
lifetime of the NLSP 

where we have inserted the value of the reduced Planck mass M P = (87tG n )" 1/2 and G N = 6.7097 • 1CT 39 GeV~ 2 
[T3 |. In order to preserve the abundances of the light elements, that are well explained by primordial nucleosynthesis, 
the lifetime of the NLSP should be shorter than tnlsp ^ 6 • 10 3 seconds [ilf . As a consequence, the latter constraint 
favours scenarios having a light gravitino, which might enter in conflict with the thermal production favouring a rather 
high reheating temperature and therefore a rather high gravitino mass, as can be seen from Eq. (|15p . 

In the case of gravitino cold dark matter, we compute the gravitino energy density fl^h 2 in our Universe as described 
above, taking into account the contributions from thermal production in the early Universe and from NLSP decay. 
In Fig. Qj] we compare the obtained gravitino relic density to the 2a range of the cold dark matter relic density of 
Eq. (|14p as a function of the gravitino mass m,Q and the reheating temperature Tr for our benchmark scenarios E to 
J. For each of the six scenarios, we also indicate the upper limit on the gravitino mass coming from the constraint 
on the NLSP lifetime, the limit between warm (WDM) and cold dark matter (CDM), as well as the region where 
the gravitino would become heavier than the NLSP. Concerning the reheating temperature, we indicate the favoured 
region with respect to leptogenesis above Tr ~ 10 9 GeV. 

Note that the contribution to il^h 2 from NLSP decay is only relevant for our point E with its neutralino NLSP 
and a rather important neutralino energy density ^nlsp^ 2 = 0.1275. The fact that this value lies already within the 
interval favoured by WMAP opens an allowed band around uIq m^o = 95.4 GeV, as can be seen in the first panel 
of Fig. [T3J For the other points, the annihilation cross section of the charged stau NLSP is more important, so that 
the resulting relic NLSP density is quite low (^nLsp^ 2 ~ 0.003 — 0.012) and its values lie below the lower limit 0.094 
of the WMAP 2<t range. 

From the graphs in Fig. [14] it becomes clear that for the chosen "collider-friendly" benchmark points, we cannot 
fulfill all three cosmological constraints at the same time. For instance, if we want a scenario featuring leptogenesis, i.e. 
having Tr > 10 9 GeV, the lifetime of the next-to-lightest SUSY particle (NLSP) would be too long for not spoiling the 
light element abundances. We therefore relax the less stringent constraint, which is the one coming from leptogenesis. 

If we then impose the constraint due to the lifetime of the NLSP, our six benchmark scenarios all lead to an 
upper limit on the gravitino mass of the order of < 10 _1 — 1 GeV. For simplicity, we propose the same value 
nig = 10 _1 GeV for all points. This respects the limit due to the NLSP lifetime, allows for gravitino cold dark matter 
with relic gravitino density that agrees with current WMAP data, and this in combination with relatively high values 
of Tr. ~ 10 7 GeV for the reheating temperature. 



V. SUPERSYMMETRIC PARTICLE PRODUCTION AT THE LHC 



In this Section, we present numerical predictions for the production cross sections of squark-antisquark pairs, squark 
pairs, the associated production of squarks and gauginos, and gaugino pairs in NMFV SUSY at the LHC, i.e. for 
pp-collisions at ^/S = 14 TeV centre-of-momentum energy. Total unpolarized hadronic cross sections 



1/2 In r 
4m 2 /S J-l/2\n 



dr / dy I ~dtf a/A (x a ,M2)f h/B (x b ,M2)^ (18) 
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FIG. 16: Examples of cross sections for charged squark-squark pair production, neutral and charged squark-antisquark pair 
production, associated production of squarks with charginos and neutralinos, and gaugino-pair production at the LHC in our 
benchmark scenario E with flavour violation in the left-left chiral sector (Arr = 0). 



are obtained through convolving the relevant partonic cross sections da/dt with universal parton densities f a /A and 
fb/B °f partons a, b in the hadrons A, B, which depend on the longitudinal momentum fractions of the two partons 
x a .b = y/Te ±y and on the unphysical factorization scales M a For consistency with our leading order (LO) QCD 
calculation in the collinear approximation, where all quark masses but the top mass are neglected with respect to 
the centre-of-momentum energy y/S, we employ the LO set of the latest CTEQ6 global parton density fit , which 
includes rif = 5 "light" quark flavours and the gluon, but no top-quark density. Whenever it occurs, i.e. for gluon 
initial states and gluon or gluino exchanges, the strong coupling constant a s (^,R) is calculated with the corresponding 
LO value of A^q - = 165 MeV. We identify the renormalization scale fiR with the factorization scales M a = Mb and 
set the scales to the average mass of the produced SUSY particles m. 

Analytic expressions for the relevant partonic cross sections can be found in Ref. [HJ . In particular, neutral squark- 
antisquark pair production can proceed from a neutral quark-antiquark pair in the initial state. At the tree-level, 
there are electroweak (s-channel photon- or Z- and i-channel neutralino- or chargino-exchange) and strong (s-channel 
gluon- and i-channel gluino-exchange) contributions. The contributing Feynman diagrams are shown in Fig. QjjJ 
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FIG. 17: Same as Fig. [16] for our benchmark scenario E with flavour violation in the left-left and right-right chiral sectors 
(Arr = All). 



where the chargino contribution in the third diagram is missing in the expressions given in Ref. [l2(]. The chargino 
exchanges, which are numerically unimportant for the dominating channels, have to be taken into account for up-type 
(down-type) quarks in the initial state and down-type (up-type) squarks in the final state. The corrected expression 
for the quark- initiated differential partonic cross section (Eq. (28) in Ref. [IH) reads 
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FIG. 18: Same as Fig. [T6]for our benchmark scenario F. 
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FIG. 19: Same as Fig. [T7]for our benchmark scenario F. 

Rather compact expressions for the appearing form factors y, Z m , Q, \yZ\ , [Gy] , \GZ\ , [GG] , [-AO*] i 

[AfZ] , [AfG] , Gnm, and Af„m (n,m = 1,2 and k, I = 1, ... ,4) and the unchanged gluon-initiated contribution to 

the cross section are given in Ref. [Hj]. The additional form factors [Cy] , [CZ\ k , [CG] k , and C™ m related to the 
chargino exchanges read 
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FIG. 20: Same as Fig. I16I for our benchmark scenario G. 



[cz] k m = 

[CGt = 



Ylx 2 w (1 - x w )s 2 
4iraa s Sjj S qql Rc 
9xw s 2 



Rc 



c tixi c T 3 i'xi y L ^ z + R ^ z > 



nm ( 1 2 2 \ 

Cqq'Z \ Ut - m q i m q j ) > 



(21) 



with all other variables defined as in Ref. [12j. We take the opportunity to also correct a few minor typographical 
errors in some of the form factors for the pair production of two up- or down- type squarks (Eq. (36) in Ref. [l2|]) 
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FIG. 21: Same as Fig. 1171 for our benchmark scenario G. 
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FIG. 22: Same as Fig. [16] for our benchmark scenario H. 
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In Figs. [TB] - 1271 we show examples of the obtained numerical cross sections for charged squark-squark pair pro- 
duction, neutral and charged squark-antisquark pair production, associated production of squarks with charginos and 
neutralinos, and gaugino-pair production at the LHC for our benchmark points E, F, G, H, I, and J and for both of 
the two considered implementations of non-minimal flavour violation in the GMSB model discussed in Sec. [Til We 
recall that the first is based on mixing between matter and fundamental messengers, leading to flavour mixing only 
in the left-left chiral squark sector and implemented at the electroweak scale through the parameter All, while Arr 
is set to zero. The second scenario involves mixing with antisymmetric messengers, giving rise to flavour violation in 
both the left-left and right-right chiral squark sectors governed by the parameter All = Arr. For the sake of better 
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FIG. 23: Same as Fig. [17] for our benchmark scenario H. 



readability, we show only the numerically most important curves as well as a selection of those that involve visible 
flavour-violating effects. 

The magnitudes of the cross sections vary from the barely visible level of 10~ 2 fb for weak production of heavy 
final states over the semi-strong production of average squarks and gauginos and quark-gluon initial states to large 
cross sections of 10 2 to 10 3 fb for the strong production of diagonal squark-squark and squark-antisquark pairs or 
weak production of very light gaugino pairs. Unfortunately, the processes whose cross sections are largest are mostly 
insensitive to the parameter All in both flavour violation scenarios, as the strong gauge interaction is insensitive to 
quark flavours and gaugino pair production cross sections are summed over exchanged squark flavours. 

Some of the subleading, non-diagonal cross sections show, however, sharp transitions in particular squark production 
channels. These transitions are directly related to the "avoided crossings" of the mass eigenvalues discussed in Sec. 
IIIII At the point, where two levels should cross, the involved squarks change character and are subject to an exchange 
of their flavour contents. Rather than the mass dependence on All, these exchanges then lead, together with the 
different parton densities in the proton, to more or less sharp transitions in the production cross sections, where the 
corresponding squarks are involved. This phenomenon is analogously observed in the case of squark and gaugino 
hadroproduction in minimal supergravity 
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FIG. 24: Same as Fig. [16] for our benchmark scenario I. 



As an example, let us discuss in detail the production of squarks and gauginos for our benchmark point E. The 
cross sections in our flavour violation scenario based on fundamental messengers are shown in Fig. [16] "Avoided 
crossings" of mass eigenvalues occur here, e.g., for down-type squarks at a value of All ~ 0.145 between the squarks 
ds and di, see also Fig. [S] Before this point, d% is characterized by a dominant sdown content, while d^ has first 
a dominant sbottom and then sstrange content. For All ^ 0.145, these contents are exchanged, i.e. c?3 is then a 
strange-squark and d^ becomes sdown-like. As a consequence, the cross sections involving the two mass eigenstates 
exchange their values, since the production of first generation squarks is preferred due to the more important parton 
density of up- and down-type quarks in the proton. This can be seen in our example for the production of down- 
type squark-squark and squark-antisquark pairs, mixed up- and down-type squark-squark and squark-antisquark pair 
production, as well as for the associated production of down- type squarks and charginos or neutralinos. For up- type 
squarks, the level-reordering phenomenon occurs at values of All — 0.09 in the range excluded by BR(6 — > sj) (left of 
the vertical dashed/green line) and is therefore not shown here. However, another effect becomes visible in the case of 
production cross sections that involve final states with up-type squarks. Some of the mass eigenstates do not present 
sharp transitions, but rather a continuous change in their flavour content. This is, e.g., the case for the lightest mass 
eigenstate u\. The corresponding production cross sections increase smoothly with the flavour violation parameter 
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FIG. 25: Same as Fig. I17l for our benchmark scenario I. 



All ) which is explained by the fact that for lower values of All the lightest up- type squark ui is mostly stop-like, but 
receives sizable contributions of the light flavours for higher All- Together with the more important parton densities, 
this results in an increase of the corresponding production cross sections. In the same way, we also observe cross 
sections that decrease with All, due to a decrease of their light flavour content. 

The same phenomena are observed in the case of our second flavour violation scenario with antisymmetric messen- 
gers, see Fig. [T7] for the benchmark point E. Note that here also "avoided crossings" between up-type squark mass 
eigenstates are observed, e.g. between U4 and U5 at All ~ 0.11, which lies, however, already in the range excluded by 
BR(6 — * S7). In this example, the M5 loses its important up-squark content to the scharm-dominatcd U4. The latter 
becomes then purely sup-like, enhancing its production cross section due to the parton density in the proton, while 
the cross sections involving M5 become less important. 

For the benchmark points H with fundamental (Fig. |2"2")) and I (Fig. |2"5")) and J (Fig. |2"T)) with antisymmetric 
messengers, we observe a third effect at All = 0.158, 0.132 and 0.114, respectively. Here, the pair production of up- 
type squark pairs (U3U4 for fundamental and U4H5 for antisymmetric messengers) exhibits an interesting resonance-like 
behaviour. It is generated by the fact that these squark mass eigenstates exchange their up and charm flavour contents 
(and also their chiralities in the case of antisymmetric messengers) at the critical ALL-values in a rather smooth way, 
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FIG. 26: Same as Fig. I16I for our benchmark scenario J. 



so that both squark mass eigenstates receive significant up- (valence-) quark contributions to their production cross 
sections in the vicinity. 

We remind the reader that in the case of flavour mixing only in the left-left chiral sector, "avoided crossings" occur 
among the qi,2, 93,4, and mass eigenstates, whereas in the case of flavour mixing in both the left-left and right-right 
chiral squark sectors, we rather observe the mass flips among the (72,3 and (74,5 mass eigenstates, respectively. Note 
also that the difference between the two flavour violation scenarios is invisible for the gaugino pair production in the 
bottom right panels of Figs. [H] - respectively, that are practically insensitive to flavour violation in the squark 
sector. 

Concerning the production of gravitinos, the cross sections achieve sizable orders of magnitude only in the case of 
a rather light gravitino, see e.g. Ref. [l4|. If the latter is too heavy, its couplings are too small to yield discoverable 
cross sections, since they are proportional to the inverse of the gravitino mass squared. In particular, this is the case 
for our scenarios with gravitino cold dark matter, where we have found a value of the order of m,Q ~ 10" 1 GeV 
derived from the different cosmological constraints. Note that in order to have a very light gravitino and consequently 
sizable production cross sections, one could consider a GMSB scenario with gravitino hot dark matter {tuq < 1 keV) 
and additional cold dark matter from stable messenger particles [43|, 0, EE EH • In scenarios with a mixing between 
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FIG. 27: Same as Fig. I17I for our benchmark scenario J. 



messenger and matter fields, however, the stability of the lightest messenger might be lost. We therefore do not 
consider the production of light gravitinos in our GMSB scenarios with additional flavour violation in the squark 
sector. 



VI. CONCLUSIONS 



While SUSY-breaking mediated by gauge interactions may be in principle attractive, since the gauge interactions 
do not induce flavour violation as do the gravitational interactions, the parameter space of minimal GMSB models is 
today severely constrained by low-energy, electroweak precision, and cosmological constraints, in particular from the 
flavour-changing neutral-current decay b — > S7. 

We have discussed several possibilities how flavour violation may still be induced in GMSB models. Focusing on 
messenger-matter mixing scenarios with fundamental or antisymmetric messengers and flavour violation in the left- left 
only or left-left and right-right chiral squark sectors, we have established collider-friendly regions of parameter space 
that are at the same time cosmologically viable and allow for the definition of benchmark points with neutralino or 
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stau NLSPs. 

Depending on the strength of the flavour- violating parameter All > we showed that splittings and avoided crossings 
appear in the mass spectra, inducing at the same time smooth or sharp transitions in the squark flavour contents. This 
induces interesting phenomenological consequences, in particular for squark and gaugino production cross sections at 
the LHC induced by valence and sea quark parton densities that differ largely in magnitude. 

We had hoped to also be able to investigate gravitino production and decay with flavour violation at the LHC. 
This was, however, not possible, since cosmological constraints require either a relatively heavy gravitino as a CDM 
candidate, with masses of 1 -4 — 10 -1 GeV and consequently unobservably small cross sections, or additional messenger 
dark matter, which would have been incompatible with our messenger-matter mixing scenario. 
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